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ABSTRACT 

I discuss the concept of quasi-state decompositions for ground states and equilibrium states of 
quantum spin systems. Some recent results on the ground states of a class of one-dimensional 
quantum spin models are summarized and new work in progress is presented. I also outline some 
challenging open problems and conjectures. 

INTRODUCTION 

The aim of this contribution is threefold. First I would like to review the notion of quasi-state 
decomposition as it was developed and used in (Aizenman and Nachtergaele, 1993a). I then show 
how interesting quasi-state decompositions for the ground states of a class quantum spin models 
can be obtained using a Poisson integral representation of the Gibbs kernel e~@ H . Secondly I will 
try to give a short review of the results we have recently obtained using quasi-state decompostions. 
These results mainly concern a class of one-dimensional quantum spin Hamiltonians introduced by 
Affleck (Affleck, 1985). The third and last part is devoted to a discussion of work in progress and 
some open problems in the subject as well as some stimulating conjectures. 

2. QUASI-STATE DECOMPOSITIONS AND POISSON INTEGRAL 
REPRESENTATIONS 

A quantum system is determined by a C*-algebra of observables A and a one-parameter group of 
automorphisms of A representing the dynamics. Here we are mainly interested in quantum spin 
systems defined on a finite lattice T for which A = Ar = ^grA 1 ) with A{ x j = M(d, C) (the 
complex d x d matrices), for all ieT. Often T is a subset of an infinite lattice (e.g. iy) and one 
is then interested in properties of the system as T tends to the full lattice (the thermodynamic 
limit). The dynamics is usually defined in terms of a Hamiltionian H = H* G A and we will give 
several examples of interesting Hamiltonians below. A state of the system is a normalized positive 
linear functional of A. A class of physically interesting states are the equilibrium states: for any 
[3 > one defines a state up by 

fA . TrAe-P" 
MA) = Tre _ m (2.1) 

for all A £ A. In this paper by ground state we always mean a state which is a limiting point of 
(2.1) as [3 — > oo, possibly under specified boundary conditions. In most examples discussed here 
(2.1) will converge to a unique limit for finite systems and the thermodynamic limit will be taken 
afterwards. 

It has been remarked many times that analyzing the structure of equilibrium states and ground 
states of a quantum system poses extra difficulties as compared to the situation in classical equilib- 
rium statistical mechanics where the states are probability measures. Sure enough these probability 
measures can be highly non-trivial as well, but the fact that there is an underlying configuration 
space which can be visualized in a concrete way helps a lot in understanding their behaviour. 
The aim of this section is to explain how at least for some quantum spin models one can obtain a 
pictorial representation of the ground states and equilibrium states which resembles somewhat the 
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situation for classical systems. More precisely we will derive a decomposition of up of the form 

up(A) = J fi(du)E w (A) (2.2) 

where fx(du) is a probability measure on a space of configurations u, and for each u, E u is a linear 
functional of Ar- Of course such an integral decomposition will only be useful if the structure of the 
E u is considerably simpler than up itself and if we are able to analyse the measure /j,(du). Below 
we show that for a certain class of Hamiltonians there is a quasi-state decomposition to which 
one can give a stochastic geometric interpretation. The resulting stochastic geometric model has 
proved to be rather helpful in revealing structural properties of the states (see Section 3 and 4). 
The linear functionals E u will in general not be states and for infinite systems they are in general 
only defined on a dense subalgebra of the form 

Aoc = \J Aa (2.3) 

Acr 

where the union is over finite subsets A of T. Implicit in (2.3) is the natural embedding of ,4a in 
Aa> for A C A'. Although the E w are typically not positive and not even bounded, they will be 
states of an abelian subalgebra of A. We therefore call these functionals quasi- states. 



2.1. Poisson integral representations 

In this subsection we consider a quantum spin model defined on a finite set of sites T, which is 
completely arbitrary at this point. The set of "bonds", denoted by B, is a collection of subsets of 
r. Soon we will specialize to the case where the bonds b G B are pairs of sites and in the case of 
antifcrromagnetic models it will be important that T be given a bipartite struture: r = Ta U Tb, 
Ta n Tb = 0- ^a and Tb are called the A- and the £>-sublattice respectively. In that situation 
a bond b = {x, y} G B will be called ferromagnetic if x and y belong to the same sublattice and 
antiferromagnetic if x and y belong to different sublattices. We denote the set of ferromagnetic 
bonds by Bf and the set of antiferromagnetic bonds by Baf- 

A complex Hilbert space H x is associated with each site x G T. Although this is not necessary 
for most of what follows, we will assume that all TL X have the same finite dimension d: TL X = <C d . 

A spin S system has d — 25 + 1, S = 5, 1, |, The observables at site x are the elements 

of B(H X ) = M(d, C). The Hilbertspace of the system is then Tir = ® xe r ^ an( ^ * ne algebra 
of observables is Ar = (SWr M(d, C). Physically important observables are usually expressed 
in terms of the spin matrices S 1 ^ 2 , S 3 which are the generators of the d-dimensional irreducible 
unitary representation of SU(2). They satisfy the commutation relations 

[S a ,SP}=iJ2^S 7 (2-4) 

7 

where a, (3, 7 <E {1,2,3} and e a/ 3 7 is the completely antisymmetric tensor with £123 = 1. One 
often combines the three spin matrices into a vector S = (S 1 , S 2 , S 3 ). The magnitude of this 
"spin- vector" is related to S by S • S = S(S + 1)1. 

A quantum spin modeled is defined by its Hamiltonian H = H* e ^4r which we write in the form 

H =-^J fe (/i 6 -l) (2.5) 

beB 

with hb G Ab = ® Ig (, M(d, C), and we assume Jb > for all beB. 

To obtain a quasi-state decomposition of up, we start from the following Poisson integral formula: 

e~ 0H = J p J p (du)K(u) (2.6) 

where: 

• pi(du) is the probability measure of a product of independent Poisson processes, one for each 
bond b G £>, running over the time interval [0,/3], and with rates Jb- For the time being we draw 
the configurations u for this process as in Figure 1 (for evident reasons the figures are drawn for 
the simplest case of a one-dimensional system: T = [a, b] G Z and B = {{x, x + l}\a<x<b — 1}. 
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• K(uj) is a product of operators hb, one for each bond occurring in u> and ordered according to 
the times at which they occur: K(u>) = Y\*hb 1 hb 2 ■ ■ ■ hb n if uj is the set of time-indexed bonds 
{(b 1 ,t 1 ),(b 2 ,t 2 ),...,(b n ,t n )} with ti < ti < ■■■ < t n . 

An important quantity is the partition function Zp — Tre~^ ff , which by (2.6) and linearity of the 
trace is given by: 

Z = [ pj(du) TtK(cu) (2.7) 



We have found that for a quite large class of interactions h the following is true: 

• TrK(ui) > for all w and this number can be computed in terms of relatively simple geometric 
properties of w. 

• the diagonal matrix elements of the operators K(uS) are all non-negative in a certain tensor 
product basis of the Hilbert space of the system. 

Let us consider some elementary examples of this before we proceed. For a more detailed discussion 
of these examples we refer to (Aizenman and Nachtergaele, 1993a). 



Figure 1. A typical configuration of the Poisson process pp(dcu). 



Example 1 

Let h be the operator which interchanges the states of the two sites, i.e. hip ® -0 = ip ® <P f° r 
any two vectors ip, tp £ <D d . In any basis of <C d the matrix elements are or 1 and a fortiori non- 
negative. K(u>) represents a permutation tt(uj) of the sites in T, and its trace is easily seen to be 
d# cycles in „ K 

The number of cycles in u> becomes a geometric property of the configuration if 
we replace the Poisson- "beeps" by two horizontal lines that cross each other as in Figure 2. With 
the convention of periodic boundary conditions in the vertical direction the number of cycles in 
the permutation tt(uj) is then equal to the number of loops in w, which we will frequently denote 
by l(cu). 

With d = 2 this interaction is equivalent to the usual spin 1/2 Heisenberg ferromagnet: h{ x y ^ = 
\ + \<*x ■ 0y, where a is the vector with as its components the three canonical Pauli-matrices. 
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Figure 2. A typical configuration of the Poisson process Pp{dui) decorated for the ferromagnetic 
models of Example 1 . 



Example 2 

For the second example (which includes the spin 1/2 Hcisenberg antiferromagnet) we require Y to 
be bipartite and all bonds are assumed to be antiferromagnetic in the sense of above (so, B = Baf)- 
We consider the interaction h defined by 

h = (-l) m ~ m ' \m, -to) (to', -to'| (2.8) 

m.m' 

where {|to}} is a basis of <C d given by the eigenvectors of the third component S 3 of the spin with 
eigenvalues to. In the case d — 2 this is the interaction for the standard Hcisenberg antiferro- 
magnet (h = | — \o x ■ a y )but the same expressions define an interesting interaction for any finite 
dimension (any magnitude of the spin) and was first studied by Affleck (Affleck, 1985) and later in 
(Affleck, 1990; Klumpcr, 1990; Batchelor and Barber, 1990). The interaction is proportional to the 
projection operator onto the singlet vector for a pair of spins and hence antiferromagnetic. It can 
be shown that the diagonal matrix elements of the corresponding K(uj) (but not the off-diagonal 
ones) are non- negative and there is again a simple formula for the trace. Each Poisson- "beep" 
is now replaced by two parallel horizontal lines as shown in Figure 3. Again this turns u> into 
a configuration of loops (assuming periodic boundary conditions in the vertical direction) and 
TtK(lu) = d*(<"). 

The examples given above are not the most general ones that can be treated but they are in some 
sense the two basic ones as will become clear at the end of this section. Unlike their classical 
analogues, quantum ferro- and antiferromagnets behave in a very different way and it is therefore 
not a surprise that they lead to two very different stochastic geometric models. 
We complete the quasi-state decomposition of by establishing the relation between expectation 
values of observables for the quantum spin system on the one hand and probabilities of events 
(or more generally expectations of random variables) in the probability measure describing the 
stochastic geometric model on the other hand. From (2.6) it follows that for any local observable 
A for the quantum spin chain 

Tr Ae~@ H f 

MA) = Tie _ fm = J »(du>)E u (A) (2.9) 

where 

f i(du>) = {Z p )- 1 p J p(du)TrK(w) (2.10) 
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and 

e ^ = ^kU (2 ' n) 

n{dto) is a probability measure on the configurations u> and for A fixed is a random variable. 

We found that for many important observables A this random variable can in fact be given a simple 
geometric interpretation. Take e.g. A — S^Sy. One then finds: 

TP ( <? 3 Q 3 \ — / C(S)l[(x, 0) and (y, 0) are on the same loop ] for Example 1 , , 

AA*x<V- j (-l^-^C^I^O) and (y,0) are on the same loop] for Example 2 1 ' 

with C(S') = (Ern=-s m2 )/( 2S ' + 1 ) = S , (S' + l)/3 and where (x,t) G Tx [0,/3] denotes a space-time 
point. I[ • ] denotes the indicator function of the event described between the brackets. Hence, 
the spin-spin correlation is proportional to the probability, with respect to the effective probability 
measure fi(dco) on the space of loop configurations, that two sites are on the same loop of ui: 

, „ 3 „ 3 , _ J C(5')Prob A1 ((a;, 0) and (y, 0) arc on the same loop) for Example 1 ,„ n , 

( 1 »' _ \ (-l) lx - yl C(S)PTob f ,((x,0) and (y,0) arc on the same loop) for Example 2 ^ > 

For u> fixed (2. If) defines a functional of the algebra of observables which is a state on a certain 
abelian subalgebra (see (2.20)). We therefore call (2.9) a quasi-state decomposition of the state 

UJp. 



Figure 3. A typical configuration of the Poisson process Pp(dw) decorated for the antiferromagnetic 
models of Example 2. 

Examples f and 2 can be combined in order to obtain a quasi-state decomposition for a class of 
models with mixed ferro- and antiferromagnetic interactions. 

In the statement of the following theorem some extra notation is used: for any configuration 
v = (< J x)xer we denote by \a) the vector in Hr given by 

k)=®kx> (2.14) 
xer 

where \a x ) is the eigenvector of S% belonging to the eigenvalue a x E {—S, —S + 1 , . . . , S}. a(t) = 
{& x (t))xer is a piecewise constant function on T x [0, j3] taking values in {—5, —S + 1, . . . , S}. A 
configuration of the Poisson process becomes a configuration of loops if we draw the ferromagnetic 
bonds as in Example I (Figure 2) and the antiferomagnctic ones as in Example 2 (Figure 3). As 
before, l(u>) denotes the number of loops in u>. We call a(t) consistent with a configuration of loops 
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u) if the discontinuities only occur at the time-indexed bonds in lo (but u> need not be discontinuous 
at a bond) and if in addition for each time-indexed bond b — ({x, y}, t) E B x [0, (3): 

, . J <j y (t + ) if {x, y} is a ferromagnetic bond 
^ [ — (T y (t~) if {x, y} is an antiferromagnetic bond 

, + s _ J &y{t~) if {a;,y} is a ferromagnetic bond 1 
x ^ \ — a y (t + ) if {x, y} is an antiferromagnetic bond 

Theorem 2.1. For a Hamiltonian H of the form (2.5) and such that hb is the interaction of 
Example 1 if b 6 Bf and is the interaction of Example 2 if b € Baf we have the following 
relations: 

i) 

(a' | e-f> H | a) = sign(*> , a) | p£(<L;) £ I[a(/3) = a', «r(0) = a] (2.17) 

<j consistent with uj 

where sign(a',cr) = (— 1) Naf W with Naf(k) denoting the number of transpositions t x>y with 
{x, y} e in any permutation it of the sites in T such that o' x = avfrr) (the sign is well-defined 
because of the bipartite structure of the model), 
ii) the partition function is given by: 

Z = Hre-f* H = J pf m (M (25 + (2.18) 

Hi) the equilibrium expectation values of observables which are functions of the operators S x can 
be expressed as 

Tr/((^) xer )e-^ 



Jfi(du;)EM) (2-19) 



where fi(duj) = Z fj 1 p[ , ( g](da;) (25 + 1)'^' and the expectation functional E u {f) is obtained by 
averaging, with equal weights, over all the spin conhgurations consistent with w: 

^ (/)= (25+V) E /("(* = «)) ( 2 - 2 °) 

cr consistent with w 



2.2. Quasi-state decompositions for more general models 

We now show that the method yielding the quasi-state decompositions of Theorem 2.1 can be 
extended, with little effort, to a much larger class of Hamiltonians of the form 

2s 

b€B {x,y}eBk=0 

It is in fact useful to think of the bonds here as a pair {x,y} of sites together with a label 
k e {1,25}. So, for a nearest neighbour model, we now have 25 bonds associated with each 
pair of nearest neighbour sites and also 25 independent Poisson processes with rates Jk(x, y) > 0. 
The Jk(x,y) are not the same as the Jk(x,y) because it will be necessary to reorganize the 25 
interaction terms in (2.21) as follows: 

2s 2s 

Mx, y)(S x ■ Sy) k = -J2 Jk(x, y)Q k (S x ■ S„) (2.22) 

k=0 k=0 

where Qk are particular polynomials of degree k to be specified below. Hence the Jk will be related 
to the Jk by a linear transformation. 

The most straigthforward way to obtain a stochastic geometric representation for a quite extensive 
class of spin 5 Hamiltonians of the form (2.21) is by representing them as a spin 1/2 system. 
To this end we introduce the (up to a phase) unique isometry V^ 5 ) : C 2S+1 — ► (C 2 )® 2S , which 
intertwines the spin 5 representation of SU(2) with 25 copies of the spin 1/2 representation: 

V^D^ig) = (£>(V2)( ff ))®2S y (s) , ge SU(2) (2.23) 
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satisfies: V^*V^ = 1 and V^V^* = P^ where P^ is the orthogonal projection onto 
the spin 5 states in the 25- fold tensor product of the spin 1/2 states. The spin 5 states coincide 
with the states that are symmetric under permutations of the 25 factors in (C 2 )® 25 , and hence 
P( s ) can be expressed as 



2S , 2S 



P (S) ® <Pi = t^t £ ® ' ^ e C ' ( 2 - 24 ) 

i=l 



i=l 

It follows by differentation of the intertwining property (2.23) that the spin 5 operators S a ,a = 
1, 2, 3, can be expressed in terms of 25 copies of the Pauli matrices (for 5=^,5" = |(7 Q ): 

V^S a = \Y j a a V^ (2.25) 

i 

We can now define a class of Hamiltomians H for a spin 5 system on a lattice T by representing it as 
a spin 1/2 model on the extended lattice T = Tx {1, . . . , 25}. With each xeTwe associate a set of 
25 sites of the form (x, k) in V. The set of bonds B is constructed as follows. For k = 1, . . . , 25 and 
b = {x,y} e 6 we^ define 6 fc = {(x, 1), (y, 1), . . . , (x, fc), (y, fc)}. Then B = {6 fc | b G B, 1 < fc < 25}. 
The interactions /i^ are defined by 

fc 

= \\ h (x,l)Xv,i) ( 2 - 26 ) 

where i s one °f the following two operators, depending on whether & is a ferro- or anti- 

ferromagnetic bond: 

h J 5 + ' a (y,i) if & is ferromagnetic 

«(x,0,(3/,0 - 1 i _ 



— \ a (x,i) ' °"(y,0 if & is antiferromagnctic 



It is then straightforward to derive (as an extension of Theorem 2.1) a quasi-state decomposition 
and the corresponding stochastic geometric representation for a Hamiltonian of the form 

\b£Bk=l ) 

where V = <S> xer 14 ■ We will describe to resulting stochastic geometric representation in detail 
below. But first we address the question of how to express a Hamiltonian of the form (2.28) in 
terms of the spin 5 matrices. In particular what are the polynomials Qk appearing in (2.22). That 
(2.28) is equivalent with (2.22) for some particualar choice of the polynomials Qk follows directly 
form the SU(2) invariance which one can check using the intertwining relation (2.23). Explicit 
knowledge of the Qk is important if one wants to check whether some particular Hamiltonian (e.g. 
one that is discussed in the physics literature) has a stochastic geometric representation of the kind 
developed here or not. 

So, we arc looking for two sets of polynomials {Q F } and {Qk F } of degree < 25 such that 



QZ(S x -S v )=V*Y[t {XilUy , l} V (2.29) 
i=i 
k 

Q£ F (S x -S y )=V*l[p {x!lUy , l) V (2.30) 



where 



t(x,l),(y,l) = 2 + l a (x.l) ■ <T(y,l) (2-31) 
P(x,l),(y,l) = 2~ l a (x,i) ' a (y,l) (2-32) 

The sites x and y in (3.29)-(3.30) just indicate any two distinct sites and play further no role. We 
therefore omit these indices in the discussion that follows. We also denote by V and P the tensor 
products of the isometry and the projection P^ over the relevant sites. 
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The following two lemma's will be useful. 

Lemma 2.2. For any k, I, k', I' e {1, . . . 25}, the following relations hold: 

= 1 Pk,l = 2 Pk,l 

Pt k ,lt k ',l = Ptk,l Ppk,lPk',l = PPk,l 

Ptk.ity itk.v = Ptk.i Pp k ,iPk>,iPk,i> = Ppk.iPk',v 

ptk ' iP = ?2^ £ pt ^ Pvk - iP = s ppi ' j 

\ ' i,j = l '>•? 

Proof :. The set of relations for the t's follows directly from (2.24) and the fact that tk,i is the 
transposition interchanging the k th factor of the 25-fold tensor product associated with the first 
site and the I th factor of the 25-fold tensor product associated with the second site. The relations 
for the p's follow from the relations for the t's by using 

Pk,i = i-t k ,i (2.33) 
which in turn follows directly from (2.31-32). 



Lemma 2.3. The operators Q F = Q F (S ■ Si) and Q£ F = Qjf F (S • Si) dcGncd in (2.29-30) 
satisfy the recursion relations 

(25 - kfQ F k+l = {(25) 2 Qf - 2fc(25 - k)}Q F k - k 2 Q F _ x (2.34) 
(25 - fc) 2 Q#i - {(2S) 2 Qt F - 2fc(45 - k + 1)}Q AF (2.35) 

for 1 < k < 25 - 1 and 

Q F = 1 Q$ F = 1 

n F - 1 i+ 1 s s n AF - : i 1 s s (2 ' 36) 

^ 1 ~ 2 252 ^ 1 ~ 2 252 

Proof :. We derive the recursion relation for the ferromagnetic case. The relation for the 
antifcrromagnetic operators can be derived in a similar way. We use the formulas of Lemma 2.2 
and the properties of V and P. 

QiQl = v*t hl v*vt hl ---t k , k v 
1 25 

= 7^2 E V*Pti,jtiA---tk,kV (2.37) 

The (25) 2 terms in the RHS of the last equality can be taken together into three groups of terms 
which are equal: (25 — k) 2 terms where i > k and j > k, k 2 terms where i < k and j < k, and the 
rest which are 2(k — 1)(25 — k + 1) terms where either i<k<jorj<k<i. The result is: 

r>F n F _ ( 2S - k f v *. v 

Vi — (25) 2 *2j2 L tk+ijk+i v 

1 2fc(25-fc + l) 
k 2 



-V*t t- V 



(25) 

where {ii, 12, ■ ■ ■ , ik+i} and {ji, J2, • • ■ , Jfe+i} are two sets of distinct indices taken from {1, . . . , 25}. 
This proves the recursion relation for the Q F . The recursion realtion for the Q k F is derived in the 
same way with one small difference: when i,j < k the cases i = j and i 7^ j lead to different terms 
due to the relations of Lemma 2.2. 



8 



For k = one has Q F = Q$ F = V*V = 1. For k = 1 we use the relations (2.25) and (3.31-32): 

Q( = V*t 1 , 1 V = V*P-^J2t iJ V 



25 2 



(2.39) 



and = I-Qf by (2.33). 

■ 

In particular this lemma shows that Q F and Qt F are two complete independent families of rotation 

29+1 2S4-1 

invariant operators on C <g) C because the recursion relations imply that the they are 
polynomials of degree k in the Hcisenberg interaction So ■ Si. The recursion relation can also be 
used to obtain explicit expressions for the polynomials Q F and Q k F - F° r 5=1 one finds: 



Qo =i 

Qi = | + ^S • Si 

Q F = -I + S -Si + (S -Si) 



or = 1 

= 5 — ^So • Si 
0^ = -1+ (S -Si) 2 



(2.40) 



For later reference we also give here a general formula for the Qj; : 

25 



(25- A;)! 



(25)! 



n 

i=25-fc+l 



(2.41) 



where zi = + 1) — 5(5 + 1) arc the eigenvalues of So • Si. From (2.41) it immediately follows 
that Q£ F (S ■ Si) is of the form 



2S-k 



Qt F (S ■ Si) = ^iP (l) 



(2.42) 



where is the orthogonal projection onto the states of total spin I in the tensor product of two 
spin 5's. One can in fact easily determine the by computing the eigenvalue Q£ F (S - Si) takes 
on a state of total spin I: 



aid 



= 2 k 



(25 -k)\ 
(25)! 



25 



] } (Z m - Zi) 



(2.43) 



m=2S-k+l 



It follows immediately that aki — if / > 25 — k and that aki > if / < 25 — k. From this 
observation it also follows that the class of integer spin antiferromagnetic Valence Bond Solid 
(VBS) chains studied in (Arovas et al., 1988) and (Fannes et al., 1989, 1992) do not alow for a 
stochastic geometric representation of the kind discussed here; the measure fi would not be non- 
negative. On the other hand the interactions of VBS models are tuned such that exact cancelations 
in the measures occur, ie., positive and negative constributions in the measure cancel each other in 
such a way that some types of fluctuations are greatly reduced. The simplest non-trivial example 
of a VBS model is the AKLT chain, introduced in (Affleck et al., 1988). For this model one can 
derive an extremely simple quasi-state decomposition for its unique ground state which reveals 
beautifully the structure of the state (Nachtergaele, 1993b). 

We now proceed to describe the stochastic geometric representation which is the result of the 
discussion above. We consider a spin 5 model on a lattice T with Hamiltonian 



25 



25 



H = - ]T J2j b , k (x,y)Q F (S x -S v )- ^ J2j bik (x,y)Q£ F (S x -S y ) 



(2.44) 



{x,y}eB F k=0 



{x,y}eB A F fe=0 



9 



It is assumed that the lattice has a bipartite structure in the sense discussed in Section 2.1. Let 
p(du>) denote the Poisson process on B as defined in Section 2. As in Examples 1 and 2 of Section 2 
Tr K(uu) can be computed in a geometric way. One proceeds as follows. Because V is an isometry 
and the cyclicity of the trace it is obvious that 

Tr^H = Tr Wr ]Y Vh b „V 

= Tr Wf PK(u>) (2.45) 

where the h bk arc defined in (2.26), J\* denotes a time-ordered product, and 

If there would be no projection operators P in (2.46), it would immcdately follow from Theorem 
2.1 that Tr K(uS) = 2 l ( w h At each site of the original lattice T there would then be 25* vertical 
lines, one for each site in T. Each lines carries a spin-^ degree of freedom. The time indexed bonds 
b k are graphically represented by the diagrams shown in Figure 4. 



Figure 4. The diagrams representing the interactions Q k and Q k . The shaded boxes represent 
the permutations TT V of the lines in each vertical segment: the index i\ at the top of the leftmost 
line indicates that this line is identified with the line number i\ under the shaded box etc. 

The projection operators P reflect the fact that only states which arc symmetric under permuta- 
tions of the 25 factors are relevant. From (2.24) it is quite obvious how the inserted projection 
operators can be incorporated into the stochastic geometric representation. It is sufficient to av- 
erage over all permutations of the lines between any two time-indexed bonds, i.e. we introduce 
a random permutation at each vertical segment in u>. By a vertical segment v in uj we mean a 
maximal time interval at a site xeT such that no bond b k = ({x,y}, k) occurs in u> during that 
time interval. For each u> we define (d7r) to be the uniform probability measure on configurations 
( 7T v)v£uj OI permutations n v of the set of 25 lines running through the segment v. It is then obvious 
that 

TtK(uj) =J ^(d 7 r)2 i « w ' 7r » (2.47-) 

It is convenient to define p~p(dQ) — pp(duj)v LU (dTr) with Cj = (u>,n). 

Theorem 2.4. For a Hamiltonian H of the form (2.44) with all Jk(x,y) > for 1 < k < 25. 
Without loss of generality we may assume that Jq(x, y) = 0. Then: 
i) 

2S 2S 

(a' | e-P« \a)=J p J p (du) ]T sign(a\ a)I[f^ HP) = o> , ^ 5(0) = a] (2.48) 

a consistent with Hi fe=l k=l 
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where u = {(T( x ,k)), &(x,k) € {~ \i \} an d & = {cr x ), & x £ {— S, S + 1, . . . , S}, and sign(a' , <r) is the 
sign defined in Theorem 2.1 (note that this sign depends only on a and o' .) 
ii) the partition function is given by: 

Z = Tre- 0H = J pj{Auj) (2.49) 

Hi) the equilibrium expectation values of observables which are functions of the operators S% can 
be expressed as 

™y -" =/**>)&(/) (2.50) 

where /i(dw) = Z^ 1 p^(duj)2 1 ^ and the expectation functional Ec(f) is obtained by averaging, 
with equal weights, over all the spin configurations consistent with uj: 

2S 

Mf) = ^j E /(!>(* = o» ( 2 - 51 ) 

cr consistent with Q fc=l 

iv) For a general local observable A one has 

up{A) = ( p(du)E a (A) (2.52) 



with 

E {Ut „ ) (A) = 2- l «""»Tr f VAV*]V h bh *v (2.53) 

There is a special case of the models (2.44) for which the stochastic geometric picture simplifies a 
lot. They are models on a bipartite lattice with only antiferromagnetic bonds (Bf = 0) and for 
which only J2S > and all other coupling constants vanish. Then only the diagram D^f appears 
and the 25 lines belonging to one site always act in unison, and therefore can equally well be drawn 
as one line. Therefore we will refer to this special cases as the single-line models as oppposed to 
the general multi-line models. The integral over the permutation variables 7r w can be carried out 
explicitely and one finds 

; (d 7 r)2'« w ' 7r » = (25*+ (2.54) 



where l(u>) is now the number of loops in the configuration u> drawn with a single line at each site 
and with the diagram D^ F at each of the time- indexed bonds as if it were a spin-^ antiferromagnet. 
Up to trivial constants the Hamiltonians of this model is given by 

beB A F 

where denotes the orthogonal projection onto the singlet state at the bond b. In fact (25 + 
l)P fc ^ is identical to the interaction h given in Example 2 (2.8), and the quasi-state decomposition 
for these models can be derived directly — without reference to the general multi-line models — as 
was done in (Aizenman and Nachtergaele, 1993a). 



3. MAIN RESULTS FOR THE ONE-DIMENSIONAL SINGLE-LINE MODELS 

In this section we only consider one-dimensional single-line models, i.e. T = [-L, L] C Z and the 
Hamiltonian is 

L 

x=-L 

where all J x are strictly positive and we will mostly assume that J x is independent of x or alter- 
natingly takes on two different values. For the proofs of the theorems in this section and further 
details see (Aizenman and Nachtergaele, 1993a). 

One is mainly interested in properties of the ground states of (3.1) in the thermodynamic limit 
(L — > 00). The existence of some specific thermodynamic limits is guaranteed by Theorem of 
our study. 
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Theorem 3.0. For each L > the Hamiltonian Hl has a unique ground state ipL and if 
J x = Jx+2 for all x, the following two limits exist for all local observables A (finite algebraic 
combinations of the spin matrices S l x , i = 1,2,3 and x G 7L): 

/AX ,. {iI>l I Ml) .. (^l | Ai/> L ) 

(A) even = lim — — — — , (A odd = lim — - — -— — (3.2 

t-zz w>l i V'l) ^d°° d (fa i ^> 

A priori it is not obvious whether (-)even an d (-)odd represent pure phases or not, and also these 
two limits are not necessarily distinct. 
Our main results are the following: 

a) Long range order 

The interaction — p(°) is obviously antiferromagnetic in nature: it favours states where nearest 
neighbour spins are antiparallel. It is therefore natural to ask whether the ground states posses 
long-range antiferromagnetic order in the sense that 

(S^Sl) ~ (-l) x m 2 for x large (3.3) 

with ra/0. This kind of behaviour is called Neel order. Ignoring the quantum fluctuations such 
a state could be pictorially represented by 

••• tmtmtmnm it mn--- 

It turns out that Neel order does not occur in the models with Hamiltonian (3.1). Absence of 
Neel order is in fact expected to hold for any one-dimensional quantum antifcrromagnet based on 
the very general argument of (Mermin and Wagner, 1966). So far a rigorous proof of the general 
statement has not been found but for the family of models under consideration here we have the 
following theorem. 

Theorem 3.1. The infinite volume ground states described in (3.2) do not have Neel order: 

lim (SgS^evenMd) = 

The interaction — p(°) could however give rise to another type of long-range order. Although for 
a chain of length > 3 there is no state in which all nearest neighbour pairs of spins are exactly 
in the singlet state, there are two states in which half of the nearest neighbour pairs do form a 
singlet: all the spins at an even site could form a singlet with their nearest neighbour on the right 
or, alternatively, with their nearest neighbour on the left. This leads to two states of periodicity 2 
with respect to lattice translations. Pictorially these two states would look as follows: 



This kind of long-range order is called dimerization for obvious reasons. Again we are ignoring the 
quantum fluctuations and one can actually not expect that these are the true ground states. Also 
in cases where there is only a preference of the spins at even sites to form singlet states with their 
left or right neighbours, we call the state (partially) dimerized. 

The stochastic geometric representation permits us to prove the following result related to the 
AfHcck-Lieb dichotomy (Affleck and Lieb, 1986). 

Theorem 3.2. For the ground states of the translation invariant model (3.1) (J x = J), one of 
the following holds: 

i) either the translation symmetry is spontaneously broken in the infinite volume ground states 

ii) or the spin-spin correlation function decays slowly (non-exponential) with 

5><S 3 S*>| = +oo 

x 

In the first case, the symmetry breaking is manifested in the non-invariance of the pair correlation: 

(S • Si) ovcn ^ (Si • S 2 )cvcn = (So ' Si) od d 
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The stochastic geometric representation in fact gives us a detailed picture of the correlations in 
the dimcrized phase. We find the following behaviour: in the state where the spins on the even 
sites are more correlated with their neighbours to the right one finds that with probability 1 some 
spins on the left half-infinite chain (—00, x] form a singlet with some spins on the right half-infinite 
chain \x + 1, +00), for each even x. In the same state this probability is < 1 for odd x. 

b ) Decay of correlations 

Theorems 3.1 and 3.2 do not specify how fast (SqS^) converges to zero as x — > 00. The decay 
rate by itself is interesting information and in particular one would like to show that the decay 
is exponentially fast (existence of a finite correlation length) in the cases where the translation 
symmetry is broken. More generally one can consider the truncated correlation function of any 
two local observables A and B. For any local observable C we denote by C x the observable obtained 
by translating C over x and let suppC be the smallest interval [a, 6] in the chain such that C is 
localized in [a + 1, b — 1]. We have an estimate for the truncated correlation function 

\(A;B X )\ = (AB X )-(A)(B X ) (3.4) 

in terms of the truncated two-point function t(x; y) of an associated two-dimensional Potts model. 
For the translation invariant case the Potts model is at its self-dual point and for the staggered 
model it can be taken to be in the high temperature phase. The precise definition of this Potts 
model is given in (Aizenman and Nachtergaele, 1993a). The comparison theorem reads as follows. 

Theorem 3.3. Let (•) denote the expectation in the ground state of a finite chain containing 
an even number of sites, or in one of the limiting states (-) C vcn or (-) dd- Then, for any pair of local 
observables A and B of the quantum spin chain there are constants Ca and Cb such that 

\(A;B x )\<CaC b Yl t ^ 

1/GsuppA 
z£suppB;c 

In this theorem suppA is roughly equal to the set of sites in the lattice on which the observable A 
acts non-trivially. 

c ) The spectral gap 

One says that the system has a spectral gap of magnitude at least 7 > in an infinite volume 
ground state (•) if for any local observable A the following inequality holds 

lim(^[i/ i ,yl])> 7 ((^A)-|(yl)| 2 ) (3.5) 

This is equivalent to the GNS Hamiltonian having a spectral gap separating the ground state from 
the rest of the spectrum. It is therefore obvious that the existence of a spectral gap would also 
follow from an estimate of the type 

\(Ae- tH B)\ <C A C B e-T (3.6) 

One can actually prove the analogue of Theorem 3.3 for the observable B being translated in 
imaginary time instead of in space and this implies the following theorem. 

Theorem 3.4. Whenever the two-dimensional Potts model associated with the quantum spin 
chain has an exponentially decaying truncated two-point function there is a spectral gap in the 
ground states (-) cvcn and (-)odd- 

4. RESULTS, OPEN PROBLEMS, AND CONJECTURES FOR THE MULTI-LINE 
MODELS 

Also here we retrict ouselves to one-dimensional models. Some important questions related to 
higher dimensional models were briefly discussed in (Nachtergaele, 1993a). We will also only 
consider translation invariant models, i.e. the coupling constants Jk do not depend on the position 
in the chain. So, the Hamiltonians under investigation here are 

L-l 2S 

H = J2 J ^ F (^-S X+1 ) (4.1) 

x——L k=0 
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where Q^ F is the polynomial of degree k defined in (2.41) and we assume that Jfe > for 1 < k < 
25. 

4.1. The excess spin operator 

In the stochastic geometric picture of the ground state provided by the Poisson integral represen- 
tation of Section 2, it is rather clear intuitively that the quantity 

z x - 2^1 a v 

y>x 

might indeed make sense when the spin-spin correlations have sufficient decay. We call these 
quantities the excess spin operators. We would like to give meaning to this infinite sum in such a 
way that the associated random variable satisfies 

E a &)= £ E K,*) ( 4 - 2 ) 

loops in u> (y,fc)on 7 

intersecting both (-oo.xjand (x,+oo) 

A crucial property then is that the total third component of the spin (or any component for that 
matter) on the collection of sites where a given loop intersects the t = axis, is always identically 
zero. So, therefore only these loops that intersect both (— 00, x] and [x + l,+oo), appear in the 
sum (4.2). Because of (2.51) decay of the spin-spin correlation function translates directly into a 
restriction on the typical size of the loops in a configuration u. 
We define local approximants S^(s) for the excess spin as follows: 

y>x 

for e > 0, x G Z, and a = 1,2,3. The S x (e) are bounded self-adjoint operators and generate 
one-parameter unitary groups, although they are not representations of SU(2). For all g £ SU(2) 
we will denote by g the vector in R 3 such that in any strongly continuous unitary representation 
U, U(g) = cxpig • S, where S denotes the vector of generators of the representation. 
For any state (•} of the infinite spin chain there is a (up to unitary equivalence) unique Hilbcrt 
space Hgns , a vector f2 e Hgns and a representation tt of the algebra of observables A^on 7Ygns 
such that 

(A) = (SI I n(A)Sl) (4.4) 

for all A e Aj_. This is called the GNS representation of the state (Bratteli and Robinson, 1981). 
A natural dense subspace of Hgns is defined by 

H loc = {TT(A)fl I A e A loc } (4.5) 

where „4i oc is the subalgebra of A formed by all finite combinations of the operators S" . 

Theorem 4.1. Let (•} be a ground state of a Hamiltonian of the form (4.1) and denote by 
Hons its GNS Hilbcrt space. If 

£ |a; 3 {S%S*}\ < +00, or for the single-line models: £ \x(S%S%)\ < +00 (4.6) 

X X 

then, for every x G Z there is a strongly continous unitary representation U x of SU(2) defined by 
the limit 

U x {g)ip = lime^-^^T/; (4.7) 

for all tp € Hgns with the following properties: 

i) U x (g) commutes with A^oo^ 

ii) U x - 1 (g) = e^< s ^U x (g) 

Hi) The domain of the generators ofU x contains Hi oc - 

The proof of this theorem will be given in (Aizenman and Nachtergaele, 1993b). 
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Corollary 4.2. There are well-defined self-adjoint generators S x of the representation U x 
which are strong limits of the local approximants (4.3) and they have a common domain for all 
ct= 1,2,3 and x £ Z. On this domain, which contains TL\ oc , one has the relation 

SS-i = + SZ (4.8) 

For situations where the excess spin operators exists one can then apply the following theorem 
to prove that sufficient decay of correlations in a ground state of a quantum spin chain with a 
Hamiltonian of the form (4.1) and S = 5, §, §, . . ., necessarily exhibits breaking of the translation 
invariance. 

Theorem 4.3. Let lo be a pure state of the spin S chain with halhntcgral S, which is invariant 
under even translations and such that the operators S x exist and have the properties stated in 
Theorem 4.1 and Corollary 4.2, then ui is not invariant under odd translations. 

Proof :. As S x is self-adjoint the following spectral projection is well-defined: 

P x = projection onto subspace where S™ £ (4.9) 

(4.8) and the S'C/(2)-commutation relations of the excess spin operators imply that the S™ have 
spectrum C |Z and that 

P x _! = 1 - P x = P x+ i (4.10) 

and hence P x commutes with any local observable: P x £ •A!/^ x+ 2k\ ^ or au ^ £ 2. In an irreducible 
representation of the observable algebra this means that P x is a multiple of the identity and because 
it is a projection we conclude that P x = X x l with X x £ {0, 1}. By (4.10) this implies 

u{P x ) = {\ even (odd) 

v x > {0 fori odd (even) v ' 

and therefore lo is not invariant under translations over an odd distance. 



4.2. A topological index as order parameter 

The configurations of loops in the plane that occur in the Poisson integral representation of the 
quantum spin chains (4.1) have interesting topological properties which we expect to be of deter- 
mining importance for the behaviour of their ground states. We will be considering configurations 
of loops in the infinite plane, which corresponds to studying the bulk of the quantum spin chain 
in the ground state. Afterwards we also make a remark about topological effects in finite and 
half-infinite chains. 

A possible configuration of loops u> for a spin-1 chain is shown in Figure 5, and in general we consider 
configuartions which are constructed using the diagrams f , k = 0, . . . , 2S. Obviously one can 
construct configurations out of these diagrams that contain infinite lines. Such configurations are 
excluded from the considerations in this section. More precisely we impose the following condition. 

Condition 4.4. Let Co be a loop conhguration for a spin-S chain. Co satisfies the Condition 
if all loops are hnite and any point (x, t) in the plane is surrounded by only a finite number of 
loops. The set of all conhgurations satisfying this condition will be denoted by fio • The probability 
measure fi(dCo) satisfies the Condition if Prob A1 (fio) = 1- 

Here, by "surrounded by a loop" , we mean that any continuous path connecting (x, t) to infinity 
must necessarily intersect that loop. 

In a translation invariant ground state of the quantum spin chain Condition 4.4 will be satisfied if 
there is suffcient deacy of the spin-spin correlation function, e.g. Condition 4.4 is satisfied if 

^|zu;(So-S x )| <oo (4.12) 

X 
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Figure 5. A typical configuration of loops Co for a spin 1 antiferromagnetic chain. The loops are 
shown with the canonical orientation with respect to X which is used in the definition of the canonical 
winding number (4.13). 

In fact (4.12) implies that the average number of loops surrounding any given point is finite: 

5> w (sgs|)| = £ Hs 3 s s x )\ 

x>0 x<0,y>0 

> i Prob(x is connected by a loop to y) 

x<0,y>0 

> jE( connected pairs x < 0,y > 0) 

> jE( loops 7 surrounding (^,0)) 

Let Co be a configuration for which Condition 4.4 is satisfied. For any (x,t) G R 2 such that (x,t) 
does not belong to a loop in Co, we define the canonical orientation of the loops 7 G Co as the 
one for which the vertical lines located at z G Z are oriented in the positive time direction if 
sign(a; — z){— l)\- x ~ z \ = 1 ([a] denotes the integral part of a), and in the negative time direction if 
sign(x — z){— l)! 2 ^ 2 ! = —1 (see Figure 5). We then define the canonical winding number W( X;t )(d)) 
of Co with respect to (x, t) by 

W(x,t)(oj) = ^ / d0 {Xtt) (4.13) 

where 9 is defined in Figure 5. By Condition 4.4 only a finite number of terms in the sum are 
non- vanishing and of course wr x<t \(Cd) is an integer. The following crucial properties, which are 
straightforward to derive from the definition of the canonical winding number, show that the value 
of the winding number is actually a property of Co, i.e. it does not depend upon the point (x,t) 
with respect to which it is computed, except for a trival dependence on the parity of x: 

W( x , t ){u) = w {XtS) (u>) = A x {Co) (4.14) 

and 

A x+1 (Co) = 2S- A x (Co) (4.15) 

for all x G H \ Z and s,feR 

What we have found is an index Aq(Co) which takes on integer values and which characterizes 
a global topological property of Co. It follows that for any Z-ergodic probability measure [i{dCo) 
on Qq, A (Co) will take a fixed value with probability 1 with respect to fi(dCo). We claim that 
the dimcrization transition observed in certain antiferromagnetic quantum spin chains can be 
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interpreted as a transition between phases with a different value for the index Aq(lu). It can be 
shown that for each integer z there are indeed Co £ O for which A (u>) = z, with the exception of 
the case s = i, where A (£j) can only be or 1. But we expect that only relatively small values 
will occur in the ground states of the Hamiltonians (4.1) although we have not been able to prove 
this so far. A second important open problem is to show that the winding number, which under 
certain conditions on the decay of spin-spin correlations is certainly measurable with respect to 
/x(dw), is also observable for the quantum spin chain. It is indeed not trivial that there exists an 
observable A G A of the quantum spin chain of which the expectation value oj(A) would be equal 
to (or determine uniquely) the value of the topological index Ao in the probability measure /x(do)) 
related to ui by (2.52) (with (3 = +oo). 

My main interest in the problems mentioned in the previous paragraph stems from the implications 
they have for the general structure of the ground state phase diagram of antiferromagnetic quantum 
spin chains. (Affleck and Haldane, 1987) make an thorough analyis of that phase diagram on the 
basis of an identifaction of the low-energy spectrum of spin chains with certain conformal quantum 
field theories. Their (non-rigorous) arguments lead to some very interesting predictions which we 
can summarize as follows. For the parameters in the Hamiltonian taking values in the complement 
of some submanifold of less than maximal dimension one expects that: i) for the half- integer spin 
antiferromagnetic chains the ground states either exhibit one of a finite number of possible types of 
dimerization, in which case they have exponential decay and a spectral gap, or the ground state is 
unique and translation invariant with slow (non-exponential) decay of correlations and no spectral 
gap; ii) for the integer spin antiferromagnetic chains there also is a finite number of possible types 
of dimerization and there is also a translation invariant phase but one with exponential decay 
and a spectral gap (the so-called Haldane-phase). So, apart from special choices of the coupling 
constants where one is on the boundary between two or more different phases, one expects that the 
integer spin chains always have expontcntial decay and a spectral gap where for the half-integer 
spin chains a gapless translation invariant phase exists. Our Theorem 4.3 (combined with Theorem 
4.1 and Corollary 4.3) shows that for half-integer spin chains with antiferromagnetic interactions 
of the form (4.1) indeed there is no translation invariant phse with rapid decay of correlations. We 
believe that the different possible phases can be characterized by the value of the index Ao. Note 
that (4.15) has no translation invariant solutions for half integer S but that there is a translation 
invariant solution A x = S for integer S. It follows that the gapless translation invariant phase for 
half-integer spin chains must have a slowly decaying spin-spin correlation function so as to make 
(4.6) divergent. 

The existence of the excess spin operators also allows for a natural interpretation of the "exotic" 
edge states in finite and half- infinite integer spin chains. It was observed in Electron Spin Resonance 
experiments on the quasi-one dimensional compound NENP (Hagiwara et al., 1990,1992; Hagiwara 
and Katsumata, 1992) that the edges interact with a magnetic moment near to it as if they carry 
spin 1/2, although the system in all other respects seems to be modelled very well by a spin 1 
antiferromagnetic chain. It was noted in (Hagiwara et al., 1990) that if the ground state of this 
system resembles the one of the AKLT model (for which the exact ground state is given in (Affleck 
et al., 1987,1988)), the spin 1/2 nature of the edge states could be easily explained. Our proof 
of the existence of the excess spin operator is valid for a quite general class of Hamiltonians. If 
one assumes that it is correct that a unique translation invariant ground state for a spin 1 chain 
implies that the index A x is identically equal to 1 (or equal to S for higher integer spin), it follows 
that the excess spin is half-integral and that the edge states should have half-integral spin as well. 

Acknowledgement. The research discussed in this paper was carried out in collaboration with 
Michael Aizenman. 
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